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Abstract 

The paper considers the wave equation, with constant or variable coefficients in IR n , 
with odd n > 3. We study the asymptotics of the distribution fa of the random solution 
at time t € IR as t —* oo. It is assumed that the initial measure has zero mean, 
translation-invariant covariance matrices, and finite expected energy density. We also 
assume that fiQ satisfies a Rosenblatt- or Ibragimov-Linnik-type space mixing condition. 
The main result is the convergence of fa to a Gaussian measure \Xoq as t —* oo, which gives 
a Central Limit Theorem (CLT) for the wave equation. The proof for the case of constant 
coefficients is based on an analysis of long-time asymptotics of the solution in the Fourier 
representation and Bernstein's 'room-corridor' argument. The case of variable coefficients 
is treated by using a version of the scattering theory for infinite energy solutions, based 
on Vainberg's results on local energy decay. 

Key words and phrases wave equation, Cauchy problem, random initial data, mixing 
condition, Fourier transform, converges to a Gaussian measure, covariance functions and 
matrices 



Supported partly by research grants of DFG (436 RUS 113/615/0-1) and RFBR (01-01-04002). 
2 Supported partly by the Institute of Physics and Mathematics of Michoacan in Morelia, the Max-Planck 
Institute for the Mathematics in Sciences (Leipzig) and by research grant of DFG (436 RUS 113/615/0-1). 
Supported partly by research grants of RFBR (99-01-00989). 



1 Introduction 



This paper can be considered as a continuation of [4]. Here we develop a probabilistic analysis 
for the linear wave equation (WE) in IR™, with odd n > 3: 



n 



u{x, t) = Lu{x,t) = £ d i (a ij (x)d j u(x, t)) - a (x)u(x, t), 

i,j=l (1.1) 

u\t=o = u (x), u\ t =o = v (x), 
d 

where <9j = — — , x G IR n , t G IR. We assume that the coefficients of equation (1.1) are constant 

OXi 

outside a bounded region, more precisely, Q>ij(x) = 5ij for \x\ > const. Moreover, we assume that 
a nontrapping condition is satisfied, i.e. all rays of a geometrical optics associated with (1.1) go 
to infinity (see Condition E3 in Section 2.1). Denote Y(t) = (Y (t),Y\t)) = (u(-,t),ii(-,t)), 
Y = (Fq ,^ 1 ) = (u ,v ). Then (1.1) becomes 

Y(t) = AY(t), t G IR; Y{0) = Y . (1.2) 

Here we set: 

A-( 

A 



^=(5 n V (1-3) 



n 

where A = £ di(a,ij(x)dj) — a (x). We suppose that the initial date Y is a random function 

with zero mean living in a functional phase space 7i representing states with finite local energy; 
the distribution of Y is denoted by Hq. Denote by fi t , t G IR, the measure on 7i giving the 
distribution of the random solution Y(t) to problem (1.2). We assume that the initial covariance 
functions (CFs) are translation-invariant, i.e. 

Q%(x,y):=E(Y*(x)Y>(y))=qi!(x-y), x,yEJR n . i,j = 0,1. (1.4) 

Next, we assume that the initial 'mean energy density' is finite: 

e := E(\v (x)\ 2 + \Vu (x)\ 2 + \u (x)\ 2 ) = g n (0) - Ag 00 (0) + g 00 (0) < oo, x G IR n . (1.5) 

Finally, it is assumed that /io satisfies a space-mixing condition. Roughly speaking, it means 
that 

Yq(x) and Y (y) are asymptotically independent as \x — y\ — > oo. (1.6) 
Our main result establishes the convergence 

Ht^Hoo, t -> oo. (1.7) 

to a stationary measure //oo, that is Gaussian and supported in space Tt. 

Previously, results of this kind have been obtained in [12, 16, 17, 18], for translation- 
invariant initial measures /xq- However, the detailed proofs were not published because of their 
length. Another drawback was the absence of a unifying argument indicating the limits of the 
methods. In this paper, such an argument is presented, which also improves the assumptions 
and makes the proofs much shorter. Like [4], the argument is based on a systematic use of 
Fourier transform (FT) and a duality argument in Lemma 5.1 (cf. works [2, 21, 22] concerning 
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FT arguments for lattice systems). This is used in conjunction with the strong Huyghen's 
principle for the WE and the fact that the rank of the Hessian of the restricted dispersion 
relation (5.20) equals n — 1. We also found a more efficient method to estimate higher order 
momentum functions and to complete some details in the proof of scattering theory results for 
the case of variable coefficients. The mixing condition has been used in [1, 2, 21, 22] to prove 
the convergence for various classes of systems. In this paper it is used in the context of the 
WE. 

We prove relation (1.7) by using the strategy similar to [4]. At first, we prove (1.7) for the 
equations with constant coefficients a^^x) = 5ij, in three steps. 

I. The family of measures p t , t > 0, is compact in an appropriate Frechet space. 

II. The CFs converge to a limit: for i, j = 0, 1, 

Qi j (x,y) = J Y\x)Y^y)p t {dY) - Qg(x,y), t - oo. (1.8) 

III. The characteristic functionals converge to Gaussian: 

tH(V) = J exp(*(r, V))tH(dY) - exp^S^, *)}, t - oo, (1.9) 

where \I/ is an arbitrary element of a dual space and <2oo is the quadratic form with the integral 
kernel (Qg(x,y)) iij=0jl . 

Property I follows from Prokhorov's Compactness Theorem with the help of arguments 
from [24]. First, one proves a uniform bound for the mean local energy in measure pt with the 
help of the FT. The conditions of Prokhorov' Theorem the follow from Sobolev's Embedding 
Theorem. Property II is deduced from an analysis of oscillatory integrals arising in the FT. 
An important role is attributed to Proposition 4.1 which establishes useful properties of the 
CFs in the FT deduced from the mixing condition. 

On the other hand, the FT alone is not sufficient to prove property III even in the case of 
constant coefficients. The reason is that a function of infinite energy gives a singular generalised 
function in the FT, and an exact interpretation of condition (1.6) in these terms is unclear. 
We deduce property III from a representation for the solution in the coordinate space, which 
manifests a dispersion of waves. In particular, for the case n = 3 and uo(x) = 0, Kirchhoff's 
formula holds: 

u(x,t) = -^ j v {y)dS(y), x e H 3 , (1.10) 

St(x) 

where dS(y) is the Lebesgue measure on the sphere S t (x) : \y — x\ — t. Then the proof of 
(1.9) proceeds with a modification of Bernstein's 'room-corridor' method, well-known in the 
random processes theory. Namely, we divide the sphere of integration in (1.10) into 'rooms' 
Rt, 1 < k < N, of a fixed width d> 1, separated by 'corridors' C\ of a fixed width p <^ d. 
As the area \R$\ ~ d 3 =const, the number iV ~ t 2 , and (1.10) becomes 

U (x, t) ~ 5^Ji, (i.ii) 



where r\ is the integral over R\. The contribution of the 'corridors' turns out to be negligible. 
Assume for a moment that v (x) and v (y) are independent for \x — y\ > p . Then r\ are 
independent if p > p , and random variable u(x, t) is asymptotically Gaussian by the CLT. 
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So, the CLT emerges from (1.10) because of integration over the sphere \x' — x\ — t and 
the first power of t in the denominator. A similar geometrical structure of an integral over the 
sphere \x' — x\ — t emerges from Herglotz-Petrovskii's formulas in a general odd dimension 
n > 5. However, the extension of the argument based on (1.10), (1.11) is not straightforward for 
n > 5 as the Herglotz-Petrovskii's formulas contain high-order derivatives of initial functions. 

We cover all odd values n > 3 in a unified techniques by modifying the approach developed 
in [4] for the Klein-Gordon equation (KGE). However, for the KGE, the solution is an integral 
over the ball \x' — x\ < t. This fact allowed us to use for the KGE a rather different approach 
based on the analysis of an oscillatory integral where the phase function ('dispersion relation') 
has a nondegenerate Hessian. For the WE, the Hessian is degenerate, which requires additional 
constructions. Here we use the fact that the 'restricted' Hessian has a maximal rank n — 1, see 
(5.20). This leads to a weaker dispersion of waves comparing to the KGE. Newertheless, we still 
obtain the representation of the solution as a sum of weakly dependent random variables. Then 
(1.9) follows from the CLT. However, checking the Lindeberg condition for the WE requires 
some delicate calculations. Here, the deficiancy in dispersive properties is compensated by the 
reduction in dimension of the domain of integration due to the strong Huyghen's principle. 

All three steps I-III of our argument rely on the mixing condition. Simple examples show 
that the convergence to a Gaussian measure may fail when the mixing condition fails: if we 
take uq(x) = and vq(x) = ±1 with probability p± = 0.5, then u(x, t) = ±t almost sure. 

Finally, we prove the convergence in (1.7) for problem (1.1) with variable coefficients. In 
this case explicit formulas for the solution are unavailable. To prove (1.7) in this case, we use 
a version of the scattering theory for solutions of infinite global energy (this strategy is similar 
to [3]). This allows us to reduce the proof to the case of constant coefficients. Namely, we 
establish the long-time asymptotics 

U(t)Y = eU (t)Y + p(t)Y , t>0, (1.12) 

where U(t) is the dynamical group of Eqn (1.1), Uo(t) corresponds to the constant coefficients 
djk(x) = 5ij, and 6 is a 'scattering operator'. The remainder, pit), is small in local energy 
seminorms || ■ ||^, WR > 0: 

\\p(t)Y \\ R ^0, t^oo. (1.13) 

The scattering theory results are based on the Vainberg's estimates for the local energy decay; 
see [23]. 

Remark 1.1 i) Under our assumptions on initial measure po, initial date Y has an infinite 
energy. Therefore, the standard scattering theory, for the solutions of a finite energy (see, e.g., 
[13]), is not sufficient for our purposes. 

ii) The order of the operators in product QUo(t) in (1.12) differs from that in Uo(t)Q considered 
in the scattering theory of finite energy solutions. An asymptotics with the order U (t)Q would 
mean that Y(t) is close to a solution of the unperturbed equation. This is impossible for the 
solutions of infinite energy as they do not converge locally to zero, hence the perturbation 
terms in the equation are not negligible. 

The paper is organised as follows. In Section 2 we formally state our main result. Sections 
3-7 deal with the case of constant coefficients: main results are stated in Section 3, the com- 
pactness (Property I) and the convergence (1.8) are proved in Section 4, and convergence (1.9) 
in Sections 5,6. In Section 7 we check the Lindeberg condition. In Section 8 we construct the 
scattering theory, and in Section 9 establish convergence (1.7) for variable coefficients. Section 
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10 discusses Vainberg's estimates. In Appendix we collected the FT calculations. 

Acknowledgements Authors thank V.I.Arnold, A.Bensoussan, I.A.Ibragimov, H.P.McKean, 
J.Lebowitz, A.I.Shnirelman, H.Spohn, B.R.Vainberg and M.I.Vishik for fruitful discussions 
and remarks. 

2 Main results 
2.1 The notation 

Denote by D the space of real functions Cg°(IR n ). We assume that the following properties 
E1-E3 of Eqn (1.1) are satisfied: 

El a,ij(x) = 5ij + bij(x), where hj(x) G D; also ao(x) G D. 
E2 ao(x) > 0, and the hyperbolicity condition holds: 3a > 

1 n 

H(x, k) = -J2 aij{x)k i k j > a\k\ 2 , x,ke JR n . (2.1) 

2 i,j=l 

E3 A non-trapping condition [23]: for (x(0), fc(0)) G IR n x IR n with fc(0) ^ 0, 

\x(t)\ — > oo as t — > oo, (2.2) 

where (x(t), k(t)) is a solution to the Hamiltonian system 

x(t) = V k H(x(t),k(t)), k(t) = -V x H(x(t),k(t)). 

Example. E1-E3 hold for the acoustic equation with constant coefficients 

u{x, t) = Au(x, t), x G IR n . (2.3) 

For instance, E3 follows because k(t) = =>- x(t) = k(0)t + x(0). 

We assume that the initial date Y belongs to the phase space H defined below. 

Definition 2.1 H = if£ c (IR n ) ©^(E") is the Frechet space of pairs Y(x) = (u(x),v(x)) 
of real functions u(x), v(x), endowed with local energy seminorms 

\\Y\\ 2 R = J (\v(x)\ 2 + \Vu(x)\ 2 + \u(x)\ 2 )dx <oo, Vi?>0. (2.4) 

\x\<R 

Proposition 2.2 follows from [14, Thms V.3.1, V.3.2] as the speed of propagation for Eqn 
(1.1) is finite. 

Proposition 2.2 i) For any Yq G Ti there exists a unique solution Y(t) G C(IR, Ti) to Cauchy 
problem (1.2). 

ii) For any t G IR ; the operator U(t) : Y i— > Y(t) is continuous in Ti. 

We now introduce appropriate Hilbert spaces of initial data of infinite energy. Let 5 be an 
arbitrary positive number. 
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Definition 2.3 Hs is the Hilbert space of the functions Y — (u, v) E Ti. with a finite norm 




Let us choose a function ((x) E C£°(IR n ) with C(0) ^ 0. Denote by #f oc (IR n ), s E JR, the 
local Sobolev spaces, i.e. the Frechet spaces of distributions u E D'(TR n ) with finite seminorms 

\\u\\ s , R := \\A s (a^/^)\\LHR^ 

where A s v := F^ x ((k) s v(k)), (k) := \J\k\ 2 + 1, and v := Fv is the FT of a tempered distri- 
bution v. For ip E D define = ^ e tk ' x ip(x)dx. 

Definition 2.4 For s E M denote W = H^ s (JR n ) © H{ oc (WC). 

Using standard techniques of pseudodifferential operators and Sobolev's Embedding The- 
orem (see, e.g. [10]), it is possible to prove that 7i° = Ti C H,~ £ for every e > 0, and the 
embedding is compact. We denote by (•, ■) scalar product in real Hilbert space L 2 (IR ra ) or in 
L 2 (IR n ) (g) IR^ or its various extensions. 

2.2 Random solution. Convergence to equilibrium 

Let (fl, E, P) be a probability space with expectation E and B(7i) denote the Borel a-algebra 
in Ti. We assume that Yq = Yq(uj,x) in (1.2) is a measurable random function with values 
in (7i, B(7i)). In other words, (ou,x) h- > Yo^,^) is a measurable map Q x IR™ -> IR 2 with 
respect to the (completed) cr-algebras £ x i3(IR n ) and i3(IR 2 ). Then Y(t) = U(t)Y is also a 
measurable random function with values in {TC,B(H)) owing to Proposition 2.2. We denote 
by Ho(dY ) a Borel probability measure in 7i giving the distribution of the Y . Without loss of 
generality, we assume (f2, S,P) = (H, B(H), /io) and Y (u,x) = u(x) for fio(du) x <ir-almost 
all (u,x) E H x IR n . 

Definition 2.5 /i t zs a Borel probability measure in 7i which gives the distribution ofY{t): 

Ht{B) = fi (U(-t)B), MB E B(H), t G IR. (2.5) 

Our main goal is to derive the convergence of the measures fj, t as t — > oo. We establish the 
weak convergence of fx t in the Frechet spaces H~ £ with any e > 0: 

lh — ? /"oo as t — > oo, (2.6) 
where //<x> is a Borel probability measure in space 7i. By definition, this means the convergence 

J f(Y)fr(dY) - | f(Y)fi OQ (dY) as t ^ oo (2.7) 

for any bounded continuous functional /(F) in space 7i~ e . 
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Definition 2.6 The CFs of measure [i t are defined by 

Qi j (x,y) = E(Y i (x,t)Y j (y,t)y i,j = 0, 1, for almost all x, y G JR n x IT (2.8) 
if the expectations in the RHS are finite. 

We set V = D © D, and (F, = (F°, + (F 1 ,* 1 ) for F = (F°, F 1 ) G 7i, and 
\j/ — (xp ^ 1 ) g £>. For a Borel probability measure /i in the space H we denote by fi the 
characteristic functional (the Fourier transform of ji) 

/}(*) = y exp(i(F, *)) fi(dY), $GD. 

A measure /z is called Gaussian (with zero expectation) if its characteristic functional has the 
form 

= exp{-^Q(*,*)}, *GD, 

where Q is a real nonnegative quadratic form in V. A measure \i is called translation-invariant 
if 

n{T h B) = fi(B), VB G B(H), h G IT, 

where T^F(x) = F(x — h). 

2.3 Mixing condition 

Let 0(r) denote the set of all pairs of open subsets A, B C !R n at distance p(A, B) > r and 
o~{A) be the cr-algebra of the subsets in Ti generated by all linear functionals F i— > (F, ^P), where 
\l/ G £> with supp ^ C ^4. We define the Ibragimov-Linnik mixing coefficient of a probability 
measure /x on 7i by (cf [11, Dfn 17.2.2]) 

if(r) = sup sup 

(AB)eo(r) Ae<7(4Be(7(B) 
/io(5) > 

Definition 2.7 Measure /i satisfies the strong uniform Ibragimov-Linnik mixing condition if 

ip{r) -> 0, r -> co. (2.10) 
Below, we specify the rate of the decay. 

2.4 Main theorem 

We assume that measure /x satisfies the following properties S0-S3: 
SO yU has the zero expectation value, 

EYq(x) = 0, x G IR n . (2.11) 



(2.9) 
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51 The CFs of /io are translation invariant, i.e. Eqn (1.4) holds for almost all x, y G JR n . 

52 /i has a finite "mean energy density", i.e. Eqn (1.5) holds. 

53 Measure /io satisfies the strong uniform Ibragimov-Linnik mixing condition, and 

oo 

Jp = J r n -V /2 (r)dr < oo. (2.12) 
o 



Remark 2.8 (1.5) implies that /i is concentrated in H$ for all 5 > 0, since 

J \\\Y \\\ 2 S Li (dY ) =e J exp(-2<J|x|) dx < oo. (2.13) 

Let £(x) = —C n \x\ 2 ~ n be the fundamental solution of the Laplacian, i.e. A£(x) = 5(x) for 
x G M n . Define, for almost all x, y G IR n , the matrix- valued function 

Q^,y) = (Qg^y)).^ = (&-y))^ (2-14) 

where 

(jj \ _i/so 00 -£*% n €- < \ (215) 

v / 1,3=0,1 2 \ q - q q - Ag J 

According to [11, Lemma 17.2.3] (see Lemma 6.2 i) below), the derivatives d^q^ are bounded 
by mixing coefficient: Wa G Z" with \a\ < 2 — i — j (including a = 0), i,j = 0, 1, 

\d?qg(z)\ < 2e oV 1 / 2 (\z\), G IR". (2.16) 

Hence, (2.12) implies the existence of the convolution £ * q^ 1 in (2.15). Denote by Qoo(\l/,^) 
a real quadratic form in V defined by 

Qoo(*,*)= ]T | g^(x,|/)vl/ l (x)^(y)rfx^. (2.17) 
Our main result is the following theorem. 

Theorem A Let n > 3 be odd, and E1-E3, S0-S3 hold. Then 

i) The convergence in (2.6) holds for any e > 0. 

ii) The limiting measure /ioo is a Gaussian equilibrium measure on 7i. 
Hi) The limiting characteristic functional has the form 

ii 0O ^) = exp{~Q 0O (W^,W^)}, *GD, 

where W : V — > H' s is a linear continuous operator for sufficiently small 5 > 0, and the 
quadratic form is continuous in 7i' s . 
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2.5 Remarks on various mixing conditions for initial measure 

We use strong uniform Ibragimov-Linnik mixing condition for the simplicity of presentation. 
The uniform Rosenblatt mixing condition [20] also is sufficient together with a higher degree 
> 2 in the bound (1.5): there exists 5 > such that 

sup e(\v (x)\ 2+s + \Vu (x)\ 2+s + \u (x)\ 2+5 ) < oo. (1.8') 
Then (2.12) requires a modification: 



I 



r n ~ 2 a p (r)dr < oo, p = mm(-^— h, (2.12') 
o 2 + o 2 



where a(r) is the Rosenblatt mixing coefficient defined as in (2.9) but without fi(B) in the 
denominator. The statements of Theorem A and their proofs remain essentially unchanged, 
only Lemma 6.2 requires a suitable modification [11]. 

3 Equations with constant coefficients 

In Sections 3-7 we consider the Cauchy problem (1.1) with the constant coefficients, i.e. 

f il(x, t) = Au(x, t), x G IR™, , s 

\ u\t=o = u (x), u\ t =o = v (x). 

Rewrite (3.1) in the form similar to (1.2): 

Y(t)=A Y(t), teM; Y(0)=Y o . (3.2) 

Here we denote 

a=(_° o ;), (3.3) 

where A = A. Denote by Uo(t), t e IR, the dynamical group for problem (3.2), then Y(t) = 
U (t)Y . Set n t (B) = fio(U (-t)B), B e B(H), t G IR. The main result for the problem (3.2) 
is the following 

Theorem B Let n > 3 be odd, and Conditions SO— S3 hold. Then the conclusions of Theorem 
A hold with W = I , and limiting measure //oo is translation-invariant. 

Theorem B can be deduced from Propositions 3.1 and 3.2 below, by using the same argu- 
ments as in [24, Thm XII. 5. 2]. 

Proposition 3.1 The family of measures {/i t , t > 0} is weakly compact in H~ £ with any e > 0, 
and the following bounds hold: 

sup E\\U (t)Y \\ 2 R < oo, R > 0. (3.4) 

t>o 

Proposition 3.2 For any ^ G V, 

= J exp(i(Y,V)) fi t (dY) ^ expi — Q^^)}, t - oo. (3.5) 

Propositions 3.1 and 3.2 are proved in Sections 4 and 5-7, respectively. We will use repeat- 
edly the FT formulas (11.2) and (11.5) from Appendix. 
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4 Compactness of the measures family 

Here we prove Proposition 3.1 with the help of FT. 



4.1 Mixing in terms of the Fourier transform 

The next proposition reflects the mixing property in terms of the FT qtf of initial CFs q$ . 
Assumption S2 implies that qo(z) is a measurable bounded function. Therefore, it belongs to 
the Schwartz space of tempered distributions as well as its FT. 

Proposition 4.1 $(k) G L p {M n ) with 1 < p < 2, and 

J (\k\ 2 -^ + \k\- 2 )\qi j (k)\dk < C( V )e < oo, i,j = 0, 1. (4.1) 

Proof We check the bound for % — j — 1 (in all other cases the proof is similar). By the 
Bohner Theorem, q^dk is a nonnegative measure. Hence, 

j \ql\k) dk\ = J ql\k) dk = g n (0) < oo, (4.2) 

owing to S2. Similarly, (2.16) and (2.12) imply that 

j \k\- 2 \qf{k)dk\ = J \k\- 2 ql\k)dk = C J \x\ 2 ~ n ql\x)dx 

+oo 

< Cie j r<f 1/2 (r) dr < C 2 ^e < oo. (4.3) 
o 

It remains to prove that measure q^ l (k)dk is absolutely continuous with respect to the Lebesgue 
measure. Function 0(r) is nonincreasing, hence by (2.12) 

r"-y /2 (r) = (n - 1) f s n - 2 V l l 2 (r-)ds < (n - 1) f s n - 2 ^ 2 (s)ds < Cy? < oo. (4.4) 

JO JO 

Then (2.16) and (2.12) imply 

f POO roo 

\ql\z)\ 2 dz < Ae 2 / (p(\z\)dz = C(n)e 2 / r n ~ l V {r)dr < C x Tpe\ \ V l ' 2 {r)dr < oo. 



Therefore, ql\k) G L 2 (IR n ). □ 
4.2 Proof of the compactness of the family {/i t } 

We now prove bound (3.4). Proposition 3.1 then can be deduced with the help of Prokhorov's 
Theorem [24, Lemma II.3.1], in a way similar to [24, Thm XII.5.2]. Formulas (11.5), (11.2) 
and Proposition 4.1 imply 

E(u(x,t)u(y,t)) =:q°°(x-y) (4.5) 
= ^ / e-^~ y) [ 1 + C °f klt q°o°(k) + ^^(C(k) + ql°(k)) + 1 "^f 1 I ff (k)} dk, 
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where the integral converges and define a continuous function. Similar representations hold for 
all i,j = 0, 1. Therefore, we have as in (1.5), 

e t := gl\0) - Ag°°(0) + g°°( ) = ^ / (fc u (*) + |A,| 2 ^ 00 (A;) + $>(kj)dk. (4.6) 
It remains to estimate the last integral. (4.5) implies the following representation for 

= 1 + CO f klt q°o°(k) + ^^(€(k) + ti°(h)) + V (*0- (4-7) 

Similarly, formulas (11.5), (11.2) imply 

ti\k) = \kf l - C °* m \ »\k) - \k\^^(C(k) + ti\k)) + l - C °l m ql\k). (4.8) 

Therefore, (4.1) and (4.6) imply that e t < C\((p)e . Hence, taking expectation in (2.4), we get 
(3.4): 

E\\U (t)Y \\ 2 R = e t \B R \ < di^eolBn]. (4.9) 
Here B R denotes the ball {x G IR™ : |x| < R} and \B R \ is its volume. □ 

Corollary 4.2 Bound (3.4) implies the convergence of the integrals in (2.8). 

Bound (3.4) also implies, similarly to (2.13), that 

sup / \\Yf 5 n t (dY) < C 5 (<p)e < oo, 5 > 0. (4.10) 
This integral estimate implies the following corollary which we will use in Section 9. 

Corollary 4.3 i) Measures fi t , t > 0, are concentrated in Hs for any 5 > 0, and the charac- 
teristic functionals fi t are equicontinuous in the dual Hilbert space H'§: for all ^i, ^ 2 £ "H's, 

|/i t (*i)-/it(*2)|< / \exp(t(Y,^ 1 -^ 2 ))-l\ f i t (dY)< J \(Y,V 1 -V 2 )\Li t (dY) 

< J lll^lb-|||^i-^2||| , i ^(^)<C'( ( 5,^,e )|||^ 1 -^ 2 ||| / 5 , t> 0,(4.11) 

where ||| • \\\' s denotes the norm in H'$. 

ii) The quadratic forms Qt(^, VP) are equicontinuous in H' s : for all ^1, ^2 £ W5, 
I Q t (*!, *x) - Q t (* 2 , * 2 )| < C| Q t (* 1 - * 2 , *0 + Q t (* 2 , # 1 - * 2 )| 

<cj l (y, * ! - * 2 > 1 (1 (r, 1 + 1 (y, * 2 > I) M^) 

< C/ 1*11* • 1*1 - ^IWs ■ (\\\*i\\\'s + IWM's) fh(dY) < C^eo)!*! - * 2 |L i > 0.(4.12) 

iii) Therefore, the quadratic form Qoo(^, *) is continuous in H^. 
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4.3 Convergence of the covariance functions 

Here we prove the convergence of the CFs of measures /i t . This convergence is used in Section 
6. 

Lemma 4.4 The following convergence holds as t — > oo: 

i,j = 0,1. (4.13) 

Proof. (4.7) and (4.8) imply the convergence for % = j: the oscillatory terms there converge to 
zero as they are absolutely continuous and summable by Proposition 4.1. For i ^ j the proof 
is similar. □ 

5 Bernstein's argument for the wave equation 

In this and the subsequent section we develop a version of Bernstein's 'room-corridor' method. 
We use the standard integral representation for solutions, divide the domain of integration into 
'rooms' and 'corridors' and evaluate their contribution. As a result, (U (t)Y , \&) is represented 
as the sum of weakly dependent random variables. We evaluate the variances of these random 
variables which will be important in next section. 

First, we evaluate (Y(t), \f r ) in (3.5) by using the duality arguments. For t e H, introduce 
the operators U'(i), U' (t) in the Hilbert space H' s , which are adjoint to operators U(t), U (t) 
in Tis- For example, 

(Y,u^t)*) = (u (t)Y,v), *gd, Yen 5 , tern. (5.1) 

The adjoint groups admit a convenient description. Lemma 5.1 below displays that the action 
of groups U' Q {t), U'(t) coincides, respectively, with the action of U (t), U(t), up to the order of 
the components. In particular, U' {t) is a continuous group in V. 

Lemma 5.1 For * = e V 

Ufa)* = (0(-,t),0(-,t)), U'(t)V = (i>(;t)M-M (5.2) 

where (f>(x,t) is the solution of Eq.(S.l) with the initial datum {uq,Vq) = and ip(x,t) 

is the solution of Eqn (1.1) with the initial date {uq,Vq) = 

Proof Differentiating (5.1) with Y, ^ e V, we obtain 

(y,&^t)*) = (U (t)Y,V). (5.3) 

Group Uo(t) has the generator (3.3) The generator of U' Q {t) is the conjugate operator 

Hence, Eqn (5.2) holds with ^ = A ip. For the group U'(t) the proof is similar. □ 
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Denote <£>(-,£) = U^t)^. Then (5.1) means that 

(Y(t),*) = {YoM',t)), tem. (5.5) 

Remark The representation (5.5) plays a central role in the proof of Proposition 3.2. A key 
observation is that &(x,t) is supported by an 'inflated' cone of thickness r where r is the 
diameter of supp ^J. The last fact follows from the strong Huyghen's principle for group U' Q (t) 
which holds for odd n > 3. Therefore, the scalar product (Y , $(•, t)) is represented as an 
integral over the 'spherical slab' of width ps r. This replaces, for a general n > 3, the Kirchhoff 
integral (1.10) written for n = 3. 

Next we introduce 'room-corridor' partition of the space IR n . Given t > 0, choose d = d t > 1 
and p = p t > 0. Asymptotical relations between t, d t and p t are specified below. Define 
h = d + p and 

a* = jh, IP = a? +d, j E Z. (5.6) 

We call the slabs R{ = {x e IR n : a j < x n < V} 'rooms' and C\ = {x e IR n : 6 J ' < x n < a j+1 } 
'corridors'. Here ) , d is the width of a room, and p of a corridor. 

Denote by Xr the indicator of the interval [0, d] and Xc that of [d, h] so that X^ez(Xr(s _ 
i^) + Xc(s — j/i)) = 1 for (almost all) s e IR. The following decomposition holds: 

(F ,$(-,t)) = 5:((F , X ^(-,t)) + (F , X ^(-,t))), (5.7) 

where Xr := Xr(% n — jh) and := Xc(x n — jh). Consider random variables r{, c J t , where 

rl = (Y , X j M;t)), d = (Y , X m-,t)), j £ 7L. (5.8) 
Then (5.7) and (5.5) imply 

(U (t)Y ^) = J2(rl + cl). (5.9) 
iez 

The series in (5.9) is actually a finite sum. In fact, (5.4) and (11.1) imply that in the Fourier 
representation, $(k,t) = A' (k)$(k,t) and $(k,t) = Q' t (k)^f(k). Therefore, 

*(*,t) = T^J e- ik *g' t (kMk) dk. (5.10) 
This can be rewritten as a convolution 

$(•,*) =Kt**, (5.11) 

where 1Z t = F~ l Q' t . The support supp^ C B- with an r > 0. Then the convolution represen- 
tation (5.11) implies that the support of the function $ at t > is a subset of an 'inflated light 
cone' 

supp$(a;,t) C {(x,t) E !R n x IR + : -r<\x\-t<r}. (5.12) 

as TZ t (x) is supported by the light cone \x\ = t as n is odd > 3. The last fact follows from 
the general Herglotz-Petrovskii formulas (see, e.g. [6, (II. 4. 4. 11) ]) and is known as strong 
Huyghen's principle. Finally, (5.8) implies 

r t — c* — for jh + t<—r or jh — t>r. (5.13) 
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Therefore, series (5.9) becomes a sum 



t 



W^o,*> = E« + 4 ^*~r> ( 5 - 14 ) 

-Aft ft 



a,s h > 1. 



Lemma 5.2 Lei n > 1, m > 0, and SO— S3 no/d. T/ie following bounds hold for t > 1 and 
Vj: 

£|ri| 2 < C(tt) d t /i, £|c?'| 2 < C(tt) (5.15) 

Proof We discuss the first bound in (5.15) only, the second is done in a similar way. 

Step 1 Rewrite the left hand side as the integral of covariance matrices. Definition (5.8) 
and Corollary 2.14 imply by Fubini's Theorem that 

E\rj\ 2 = (xi(x n )xi(y n )Qo(x-y),$(x,t)®<S>(y,t)). (5.16) 

The following bound holds true (cf. [19, Thm XI. 19 (c)]): 

71 — 1 

sup |$(x,t)l = 0(t~—), t^oo. (5.17) 
In fact, (5.10) and (11.2) imply that $ can be written as the sum 

*0M) = tA^E / e-^l^a^lA;!)^) dfc> (5 . 18) 

R" 

where a ± (|fc|) is a matrix whose entries are linear functions in \k\ or l/|/c|. Let us prove the 
asymptotics (5.17) along each ray x = vt + Xq with |n| = 1, then it holds uniformly in x G IR™ 
owing to (5.12). In polar coordinates, we get from (5.18), 

$(nt + x ,t) = — V V r { ( e^^^-^^dkD^ik) dS(k)\ dr. (5.19) 
(2n) n ± Jo I J I 
v ' ± \k\=r ) 

This is a sum of oscillatory integrals with the phase functions <f)±{k) = kv ± \ k\. The standard 
form of the method of stationary phase is not applicable here as the set of stationary points 
{k G IR n : V0±(/c) = 0}, is a ray v = ±k/\k\, and the Hessian is degenerate everywhere. 
On the other hand, restricted to the sphere \k\ = r with a fixed r > 0, each phase function 
0± := 0±l|fc|=r = kv ±r, has two stationary points ±vr, and the Hessian is nondegenerate 
everywhere: 

rank (Hess (f>±(k)) = n - 1, \k\ = r. (5.20) 

Hence, the inner integral in (5.19) is 0(t~^ n ~ 1 ^ 2 ) according to the standard method of station- 
ary phase, [8]. At last, for the integral in r in (5.19), has the same asymptotics in t as fy(k) 
decay rapidly at infinity. 

Step ,2 According to (5.12) and (5.17), Eqn (5.16) implies that 12 

E\ri\ 2 < Ct' n+1 J xi(x n )\\q (x-y)\\dxdy, (5.21) 

SfxSf 
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where S% is an 'inflated sphere' {x G IR n : — r < \x\ — t < r} and \\qo(x — y)\\ stands for the 
norm of the 2 x 2-matrix, (qo(x — y)\. The estimate (2.16) implies then 



E\rl\ 2 <Cr n+1 J xi(x n )<P 1/2 (\x-v\)dxdy : (5.22) 
For large t, this integral can be reduced to the product of the spheres S t = {x <E JR n : \x\ — t}: 
E\rj\ 2 < CrfH-^ 1 J xi(x n ) J V? 1/2 (k - y\) dS(y) dS(x), (5.23) 



St \s t 

where dS is a Lebesgue measure on the sphere. The inner integral can be estimated by a direct 
computation owing to (2.16): 

J <p 1/2 (\x - y\) dS(y) = C{n) J* r n_ V 1/2 (r)dr < C{n)Tp, t > 0. (5.24) 
s t ° 

Therefore, (5.23) and (2.12) imply 

E\ri? < Ci(n)r 2 r n+1 J X i(x n )dS(x) < C 2 {n) d t /t. □ 

s t 

6 Convergence of the characteristic functionals 

In this section we complete the proof of Proposition 3.2. As was said, we use a version of 
the CLT developed by Ibragimov and Linnik. This gives the convergence to an equilibrium 
Gaussian measure. If Qoo(^, = 0, Proposition 3.2 is obvious. Thus, we may assume that 

Qoo(*,*)^0. (6.1) 

Choose < 5 < 1 and 

Pt-t 1 - 5 , ck~-L, t^oo. (6.2) 
hit 

Lemma 6.1 The following limit holds true: 

N t (M + (^) 1/2 )+N?(^(p t ) + ^)^0, t^oo. (6.3) 
Proof. Eqn (6.3) follows from (4.4) as (6.2) and (5.14) imply that N t ~ Int. □ 
By the triangle inequality, 

|&(*)-/ioo(*)l < \Eexp{i(U (t)Y ,^)}-Eexp{q2 t ri}\ 

+ 1 exp{-^]T £|r|| 2 } - exp{-iQoo(^, *)}| 

+\Eexp{i£ t ri} - exp{-^2 t E \ r3 t\ 2 }\ 
= J 1 + / 2 + J 3 , " (6.4) 
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N t 

where the sum J2t stands for J2 ■ We are going to show that all summands I±, I2, h tend 

j=-N t 

to zero as t — > 00. 

Step (i) Eqn (5.14) implies 

h = |£exp{^}(exp{^} - 1)| < £ £|c?| < E^N'I 2 ) 172 - ( 6 - 5 ) 
From (6.5), (5.15) and (6.3) we obtain that 

h < CN t ( Pt /tf 2 -> 0, t^oo. (6.6) 
S^ep ^ By the triangle inequality, 

/ 2 < ^E^IWI 2 -Qoo(^,^)|<i|Qi(^,^)-Qoo(^,^)| 

= / 2 l+/ 22 + / 23 , (6.7) 

where Qt is a quadratic form with the integral kernel (Q % j(x, y^j. Eqn (4.13) implies that 
I21 — > 0. As to J 22 , we first have that 

/ 22 <E^I^I- ( 6 - 8 ) 

The next lemma is a corollary of [11, Lemma 17.2.3]. 

Lemma 6.2 Let ^ be a complex random value measurable with respect to a -algebra a(A), rj 
with respect to a-algebra a(B), and dist(^4, B) > r > 0. 

i) Let {E\i\ 2 fl 2 < a, {E^ 2 ) 1 ' 2 < b. Then 

\E£ri-E£Eri\ < Cab y? 1/2 (r). 

ii) Let |£| < a, \r]\ < b almost sure. Then 

\E£r]-E£Eri\ < Cab <p(r). 

We apply Lemma 6.2 to deduce that I22 — > as t — > 00. Note that r{ = (Yq(x), x J r {x n )(TZt* 
ty)) is measurable with respect to the cr-algebra a{R ] t ). The distance between the different 
rooms R\ is greater or equal to p t according to (5.6). Then (6.8) and SI, S3 imply, together 
with Lemma 6.2 i), that 

I22 < CNf^pt), (6.9) 

which goes to as t — > 00 because of (5.15) and (6.3). Finally, it remains to check that I23 — > 0, 
i — > 00. By the Cauchy - Schwartz inequality, 



< + c(^(E t r|) 2 ) 1/2 (iV t ^N| 2 ) 1/2 . (6.10) 
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Then (5.15), (6.8) and (6.9) imply 

£(E/i) 2 < HE\rl\ 2 + 2j2 j<:l E\rir l t \ < CN t d t /t + C^ 1 ' 2 ^) < C 2 < oo. 
Now (5.15), (6.10) and (6.3) yields 

hs < C x N 2 tPt jt + C 2 N t { Pt /t) 1 / 2 - 0, t ^ oo. (6.11) 
So, all terms J 2 i, -^22, ^23 in (6.7) tend to zero. Then (6.7) implies that 

/2<^lE^| 2 -2oo(^)H0, t^oo. (6.12) 

Step (in) It remains to verify that 

h = \EeM^2 t ri} ~ exp{-^(E^) 2 }| - 0, t - 00. (6.13) 
Using Lemma 6.2, ii) we obtain: 

N t 

\Eexp{iJ2 t r J t} ~ H Eexv{ir{}\ 

-Aft 

AT* AT* 

< \E exp{ir t -JV *} exp{i E r{} — E exp{irT~ Nt }E exp{i E H}\ 

-N t +1 -N t +1 
N t N t 

+ \Eexp{ir; Nt }Eexp{i E 4} ~ H E exp{ir{}\ 

-Nt+1 -N t 
N t N t 

< C V ( Pt ) + \Eexp{t E H}~ II Eexp{iri}\. 

-Nt+1 -Nt+1 

We then apply Lemma 6.2, ii) recursively and get, according to Lemma 6.1, 

\Eexp{iJ2 t ri} ~ H Eexp{ir 3 t }\ < CN t (p(p t ) -> 0, t -> 00. (6.14) 
It remains to check that 

At 1 

in^exp{^}-exp{--E^IW| 2 }|^0, t^oo. (6.15) 

-JV t 2 

According to the standard statement of the CLT (see, e.g. [15, Thm 4.7]), it suffices to verify 
the Lindeberg condition: Ve > 

^E t ^K| 2 - °> ( 6 - 16 ) 

Here 07 = X^l r t| 2 ; an d -^5/ = EX$f, where X5 is the indicator of the event |/| > 5 2 . Note 
that (6.12) and (6.1) imply that 

O t -> Qoc(^) ^0, t -> OO. 

Hence it remains to verify that Ve > 

E t W| 2 -*0, t^oo. (6.17) 
We check (6.17) in Section 7. This will complete the proof of Proposition 3.2. □ 
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7 The Lindeberg condition 

The proof of (6.17) can be reduced to the case when for some A > we have, almost sure that 

\u Q (x)\ + \v (x)\ < A < oo, x G IR n . (7.1) 

Then the proof of (6.17) is reduced to the convergence 

J2 t E\r{\ 4 ^0, t^oc (7.2) 

by using Chebyshev's inequality. The general case can be covered by standard cutoff arguments 
taking into account that bound (5.15) for E\r 3 t \ 2 depends only on e and ip. The last fact is 
evident from (5.21)-(5.24). We deduce (7.2) from 

Theorem 7.1 Let the conditions of Theorem B hold and assume that (7.1) is fulfilled. Then 
for any ^GD there exists a constant C(^) such that 

E\r J t \ 4 < C(V)A 4 d 2 t /t 2 , t > 1. (7.3) 

Proof. Step 1 Given four points X\, x 2 , x 3 , x± G IR™, set: 

Mq 4 \xi, x 4 ) = E (Y ( Xl ) <g> ... <g> Y (x 4 )) . 

Then, similarly to (5.16), Eqns (7.1) and (5.8) imply by the Fubini Theorem that 

E\4\ 4 = (xiW) . . . xi(x n 4 )M^(x 1 , . . . , x 4 ), $(x u t) ® . . . ® $(x 4 , *)>■ (7.4) 

Let us analyse the domain of the integration (IR™) 4 in the RHS of (7.4). We partition (IR™) 4 
into three parts W 2 , W3 and W 4 : 

4 

(IR n ) 4 = M Wi, Wi = {x = (x 1 ,x 2 ,x 3 ,x 4 ) G (IT) 4 : \ Xl - x t \ = max \x x - x p \}. (7.5) 
i=2 ' p=2 ' 3 ' 4 

Furthermore, given x = (xi, x 2 , X3, x 4 ) G Wi, divide IR™ into three parts Sj , j = 1,2,3: 
IR" = Si U S 2 U 53, by two hyperplanes orthogonal to the segment and partitioning it 

into three equal segments, where x\ G S\ and X { G S3. Denote by x p , x q the two remaining 
points with p,q^ Set: Ai = {x G Wi : x p G Si,x q G S3}, Bi = {x G Wi : x p ,x q £ Si} 
and Q = {x e Wi : x p ,x q G" S3}, i = 2,3,4. Then Wi = AiU Bill d. Define the function 
m[, 4) (x), x G (IR™) 4 , in the following way: 

™(4)/^ _ / M^ 4 \x) -qofa -x p ) ®q (xi - x q ), x G Ai, ( „ , 

m ° {X) m ~ { Mi 4 \x), xeB^ Cl . (7 ' 6) 

This determines nig 4 ^) correctly for almost all quadruples x. Note that 
(xjW) • • • X j r {xl)qo{xi - x p ) <g> q (xi - x q ), $(x u t) <g> . . . <g> $(x 4 , f)) 

= (xi(Xl)xl(Xp)qo(xi - X p ), t) <g> $(x p , *)> (X J r«)X J r(^)9o(^ - X 9 ), $(^, t) <g> *)). 
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Each factor here is bounded by C(^) d t /t. Similarly to (5.15), this can be deduced from an 
expression of type (5.16) for the factors. Therefore, the proof of (7.3) reduces to the proof of 
the bound 



$(x 4 ,t))\ < C(tf)AX/i \ t>\ 



I t := \(xl(xV . . . xUxlWoHx!, . . .,x 4 ), H>( Xl ,t) ® . . 
Step 2 Similarly to (5.21), Eqn (5.17) implies, 

/i<C($)r 2n+2 J xi(xi)...xi(x2)\m i *\x 1 ,...,x 4 )\dx 1 dx 2 dx 3 dx 4 , 



(7.7) 
(7.8) 



(sir 



where is an 'inflated sphere' {x G IR n : — r < \x\ — t < r}. Let us estimate m^ 4 ** using 
Lemma 6.2, ii). 



Lemma 7.2 For each i — 2, 3, 4 and almost all x G Wi the following bound holds 

|m[, 4) (xi, . . .,a:4)| < CA 4 ip(\x l - Xi\/3). 



(7.9) 



Proof. For x G Ai we apply Lemma 6.2, ii) to IR 2 ® IR 2 -valued random variables £ = 
io(^i) ® ^o(^p) and 77 = lo(xj) ® io(^)- Then (7.1) implies the bound for almost all x e Ai 



m£\x)\ <CAV(ki-^|/3). 



(7.10) 



For x G £>j, we apply Lemma 6.2, ii) to £ = io(^i) an d ^ — ^o(^p) ® ^o(^g) ® Yo(xi)- Then SO 
implies a similar bound for almost all x G £>j, 

\m£\x)\ = |M (4) (x) -EY ( Xl ) ®E(Y (x p ) ®Y (x q ) ®Y ( Xi ))\ < CA 4 <p(\ Xl - Xi \/3), (7.11) 

and the same for almost all x G Cj. □ 
S'tep 3 It remains to prove the following bounds for each % — 2, 3, 4 (cf (5.22)): 



:= y X J r K) • • • xl^MxMlx! - x i \/3)dx 1 dx 2 dx 3 dx 4 < Cdft 

(sir 



2 + 2n-A 



(7.12) 



where X,i is an indicator of the set W^. In fact, this integral does not depend on i, hence set 
i = 2 in the integrand. Similarly to (5.23), this integral can be reduced to the product of four 
spheres S t = 5°: for large t, 

( 



V t (t)<Cr* j rt(xl)<p(\x 1 -x 2 \/3) 

(St) 2 



J xi(xl) j X 2 {x) dS(x 4 ) dS(x 3 ) 
_s t \s t J 



dS(xi) dS(x 2 ). 

(7.13) 

Now a key observation is that the inner integral in dS(x 4 ) is 0(\x 1 — x 2 \ n l ) as X 2 (x) = for 
\x 4 — Xi\ > \x± — x 2 \. This implies 



Vi(t) <Cr" J xiK) J <p(\ Xl - x 2 \/3)\ Xl - x.r 1 dS(x 2 ) dS( Xl ) J xl(x n 3 ) dS(x 3 ). (7.14) 
S t \St J s t 
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The inner integral in dS(x 2 ) can be estimated by a direct computation: similarly to (5.24), 
<p(\xi - x 2 \/S)\x 1 - x 2 \ n ~ 1 dS(x 2 ) = C(n) / r 2n ~ 3 ip(r/3) dr 



s t 

r-2t 







< d(n) sup r n ~V V V/3) / r n - 2 ip 1/2 (r/3) dr. (7.15) 

re[0,2t] JO 

The 'sup' and the last integral are bounded by (4.4) and (2.12), respectively. Therefore, (7.12) 
follows from (7.14). This completes the proof of Theorem 7.1. □ 

Proof of convergence (7.2). The estimate (7.3) implies, since d t < h ~ t/N t , 

f2 N t 

8 Scattering theory for infinite energy solutions 

As was said in Sections 1-2, we reduce the proof of Theorem A to Theorem B by using a 
special version of the scattering theory, for solutions of infinite energy. Recall, U(t), U (t) are 
the dynamical groups of the Cauchy problems (1.2), (3.2), respectively. 

Theorem 8.1 Let E1-E3 hold, and n > 3 be odd. Then there exist 5, 7 > and linear 
continuous operators 6, pit) : Hs — > Ti such that for Y Q e 7i$ 

U(t)Y = eU (t)Y + p(t)Y , t>0, (8.1) 

and for any R > there exists a constant C = C(R, 5, 7), such that for Y e Hs 

\\p(t)Y \\ R <C e-^\\\Y \ls, t>0. (8.2) 

We deduce Theorem 8.1 with the help of duality from a special version of the finite energy 
scattering theory that is developed below. Denote ||| • |||^ the norm in the Hilbert space H' s , 
dual to TL§. 

Lemma 8.2 The following bound holds true: 

W Q {t)S>\\' 5 < Ce^ll^lll^, WeH's, (8.3) 
Lemma 8.2 follows by duality from Lemma 8.3: 

Lemma 8.3 Let E1-E2 hold. Then W5 > the operator Uo(t) is continuous in Hs, and there 
exists a constant C = C{5) > such that for Y e Hs 

|||C/o(0^b|]U < C'e <5| * l lll^o|!U, telR, (8.4) 
Proof. It suffices to consider Y £ V. Then U (t)Y Q G V. Denote 

E S (Y) = J e- 25 W(\v(x)\ 2 + \Vu(x)\ 2 ) dx, Y = (u,v) G H 5 . 
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Then the derivative 

E s (U (t)Y ) = 2 J e- 2Slxl {u(x,t)u(x,t) + Vu(x,t)Vu(x,t)) dx. 

Substituting u(x, t) = Au(x,t) and integrating by parts, we obtain 

E s (U (t)Y ) = -2 J Ve- 2S ^ ■ Vu(x,t)u(x,t) dx. 

Then \Es(U (t)Y )\ < 25Es(U (t)Y ) by the Young inequality. Therefore, the Gronwall inequal- 
ity implies 

E s (U (t)Y ) < e^E s (Y ). 

In other words, 

J e- 2S W(\u(x,t)\ 2 + \Vu(x,t)\ 2 ) dx < e 2m E 5 (Y ). (8.5) 
It remains to estimate the norm |-u(-,t)|,5, where 

\u(-,t)\ 2 = J exp (-25\x\)\u(x,t)\ 2 dx. 

In fact, we have: 

t t 
H;t)\s<\u°(x)\ s + \ Ju(;r)dr\ s <\u°(x)\ 5 + \ J \u(-,r)\ s dr\. 



Using (8.5), we get 

|«(-,Oli<Ce^l|||yo|||i. □ 

Now we employ Vainberg's bounds for the local energy decay; this plays the key role in 

the proof of Theorem 8.1. We use the Sobolev space = H 1 ^^) © L 2 (B R ) with the norm 

o 

|| • ||(7j), R > 0. Recall that H~ 1 (B R ) is a completion of D R = {ip G D : suppip C -Br} in the 
Hilbert norm of the Sobolev space H~ 1 (M n ). We will use the following convenient description 
of a dual space (see, e.g. [25, Sen 1.3.4]). 

o 

Lemma 8.4 H~ l (B R ) is the dual to the Hilbert space ^(Br) with respect to the scalar product 
(;■)■ 

Corollary 8.5 The dual space to the Hilbert space 'H(r) with respect to (•, •) is 

H[ R) =H-\B R )®L 2 (B R ). (8.6) 

Note that T~L[ R ) is a subspace of H' s with any 5 G IR. 

Definition 8.6 H! denotes the space U r> qH'^ endowed with the following convergence: a 
sequence \l/ n converges to ^ in Ti! iff^R > s.t. all \l/ n G T~t'( R ), and \l/ n converge to ^ in the 
norm ofH!^ R y 
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Below, we consider the continuity of the maps from TV only in the sense of the sequential 
continuity. Vainberg's results imply the following lemma which we prove in Appendix. 

Lemma 8.7 Let E1-E3 hold, and let n > 3 be odd. Then \/R,Ro > there exist constants 
a, C(R, Rq) > andT = T(R, R ) > such that for V e H[ R) 

\\U'(t)nL H B Ro) em { B RQ ) < C(R,R )e- at \\n{ R) , t > T, (8.7) 

where || • ||'^ denotes the norm of the dual Hilbert space TV^ R y 

Now we are in position to discuss our version of the scattering theory for solutions of finite 
energy. 

Proposition 8.8 Let E1-E3 hold, and n > 3 be odd. Then there exist 5, 7 > and linear 
continuous operators W, r(t) :Ti'^Ti' s such that for ^ e TV 

U\t)V = U' {t)W^> + r(t)V, t > 0, (8.8) 

and for any R > there exists a constant Cr = C (R, 5, 7) such that 

Mtm\' s <C R e-^\\*\\[ R ), t>0. (8.9) 

Proof. We apply the standard Cook method: see, e.g., [19, Thm XI. 4]. Fix \P e Ti-'m) and 
define W^>, formally, as 

= lim U^(-t)U'(t)y = Uoi-TjU' '(7\)* + / —U' Q (-t)U' '(t)* dt (8.10) 
t^co J dt 

Ti 

with an appropriate 7\ > 0. We have to prove the convergence of the integral in the norm of 
the space TV S . First, observe that 

jU' {t)^ = A' U>(t)*, jU'm = A'U'W, 

where A' and A' are the generators to the groups U^t), U'(t), respectively. Similarly to (5.4), 
we have 

where Au = J2di(aij(x)dju) — a (x)u. Therefore, 



f t K(-t)U'(t)V = U'(-t)(A' - A' )U'(t)V. (8.12) 



(8.11) and (5.4) imply that 



A'-A- ° A ~ A " V (8.13) 



Observe that A — A = J2dibij(x)dj — a (x), where bij(x), a (x) e C™{B Ro ) with some R > 0, 
according to El. Therefore, by (8.3), we have that 

\WoH)(A' - A' )U'(t)n's < Ce St \\\(A' - A' )U'(t)n's 

<Ce St \\((A'-A' )U'(t)*)°\\ H - HBRo) <Ce 5t \\(u'my\\ m{BRo) , t>0. (8.14) 
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Then (8.12) and (8.7) imply, for t > T = T(R,R ), that 

\\\jU>(-t)U'm\\\' s < C(R,R )e St e- at \\n[ R) = C(R)e^\\n{ R) , (8.15) 

where (3 = a — 5. Choose 5 > sufficiently small: 5 < a. Then we have (3 > 0, and (8.15) 
implies 



, a 



dt 



U' (-t)U' (t)V\\\'s dt < C^mU < oo. 



T 



Therefore, the existence of the limit in (8.10) follows if we choose Ti = T. Furthermore, the 
operator W : W — > H' s is continuous, and (8.10), (8.15) imply 

- W )*|' a < C 2 (R)e^\\n{ R) , t > T. (8.16) 

Let us now choose 5 < a/2. Then 8 < (3 = a — 5 and '-f = (3 — 5 = a — 25 > 0. Finally, set 
r(t)tf := C/'(t)tf - [/*(t)Wtf , then (8.16) and (8.3) imply 

iHtmw's = WW - u^wmw's = imt)(u>(-t)u>(t) - wmt 

< C7s(i2) e «B(^(-0^(0-W)*r«<^(i2)e- 7t ||*||' (Jl) , t >T. □ 

Proof of Theorem 8.1. (8.8) implies that for Y e H and * G for any i? > 0, 

(C/(t)y o ,*) = (C/o(i)^o,^) + (>o,r(0*), t>0 (8.17) 

By Proposition 8.8, operators W R = Wl^' and r^(t) = r(t)| W ' are continuous as maps 
7~L't R \ — > 7^. Therefore, the reflexivity of the Hilbert spaces implies the existence of the adjoint 
continuous operators Or = W' R : Hs — > and p/?(t) = r#(£) : 7^ — > 7i(R) for any i? > 0. It 
remains to define (QYq)\b r = QrYq and (p(t)Y )\B R = Pn(t)Y for any R > 0. □ 

9 Convergence to equilibrium for variable coefficients 

We deduce Theorem A from next two Propositions 9.1 and 9.2 (cf. Propositions 3.1, 3.2). 
Proposition 9.1 Family of measures {p t , t e 1R} ; is weakly compact in H~ £ , Ve > 0. 
Proposition 9.2 For every ^ E V, 

jim = J exp(t(Y,y))fi t (dY) expi-^QUV,*)}, * - oo. (9.1) 

Proof of Proposition 9.1. Similarly to Proposition 3.1 , Proposition 9.1 follows from the 
bounds 

sup E\\U(t)Y \\ R < oo, R>0. (9.2) 

t>o 

Theorem 8.1 implies that 



£||l/(*)^o||* < 2E\\QU (t)Y \\ 2 R + 2E\\r(t)Y \\ 2 R < C 1 (R)E\lU (t)Y \f s + C 2 (R) e - 2 ^ E\\\Y \\ 
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Then (9.2) follows from (4.10) and (2.13). □ 
Proof of Proposition 9.2. Let ^ G 7~C'my Then Theorem 8.1 implies that 

jl t (^) = Eexp{t(U(t)Y ,^}} = Eexp{i(eU (t)Y + r(t)Y ,V)} 

= Eexp{t(&U (t)Y ,^}} + u(t), (9.3) 

where 

u(t) = E [exp{i{QU (t)Y , *>}(exp{i(r(Olo, *)} - 1)] . 
Note that z/(t) vanishes as t — > oo. In fact, Theorem 8.1 implies as above, 

KOI < £ |(r(t)F , *>| < C(i2)||*||' (fl) f?||r(*)yo||fl - 0, t ^ oo. (9.4) 
Finally, Proposition 3.2 and Corollary 4.3 imply that 

Eexp{i(QU (t)Y ,V)} = E exp{t(U (t)Y , W*)} -> exp{-^Q 0O (^, W^)}, t -> oo, (9.5) 

and (9.3)-(9.5) imply (9.1). □ 

10 Vainberg's estimates 

In this section we prove Lemma 8.7. 

Proposition 10.1 Let E1-E3 hold, and n > 3 be odd. Then \/R, R > there exist constants 
a, C(R, Ro) > and T = T(R, Ro) > such that for Y G Ti with supp Y C Br , 

Qk 

\\-^U(t)Y \\ R < C(R, R )e- at \\Y \\ Ro , t > T, k = 0, 1, . . . (10.1) 

Proof Conditions E1-E3 allow us to use Theorem X.4 from [23] which implies an asymptotic 
expansion for the solution Y(-,t) = U(t)Y , 

N dj 

F(x,t)=^^^%(x) + ^(x,t). (10.2) 

j=l 1=0 

Here Im Uj is a nonincreasing sequence, Im Uj — > — oo, < oo, and the remainder Z^(x,t) 
satisfies the bounds (10.1) when Im u^ + i < 0. It remains to prove that all terms with 
Im ujn > vanish. Assumptions El and E2 provide an a priori bound for finite energy 
solutions. Therefore, in accordance with Theorem 8 (or Lemma 10) of [23, Ch X], all increasing 
terms with Im Uj > or with Im Uj = and I > vanish. By the same reason, each amplitude 
Yji(x) with Im Uj = and I = has a finite global energy and hence is an eigenfunction of the 
operator L (see (1.1)) with the eigenvalue —Uj. Therefore, an extension of Kato's Theorem 
[7, Thm 2.1] implies the absence of a discrete spectrum inside the continuous spectrum, i.e. 
Yji(x) = if Uj 7^ 0. Amplitude with Uj = and I = vanishes since the operator L is strictly 
negative. □ 
Therefore, we get by duality the following bounds for ^ G 7~L[ R y 



Qk 

W U\t)nH-HB Ro )eLHB Ro )<C(R : R )e- at \\n{ R ): t>T, A; = 0,1,... (10.3) 
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Recall that (5.2) implies the representation U'(t)^f = (ip(-,t),i^(-,t)) where tp(x,t) is a 
solution to = Atp. Then (10.3) with k — implies 



U(;t)\\ LHBRi)) < C(R,R )e- at \\n' (R) , t > T. 
Similarly, (10.3) with k — 1 implies that 

U(;t)\\ LH B Ro) < C(R,R )e-°*\\n{ R) , t > T. 
Also, by virtue of = Aip, we get 

\\A^(.,t)\\ H - HBRo) < C(R,R )e- at \\n{ R) , t > T. 



(10.4) 



10.5) 



(10.6) 



Note that (10.5) is a part of the bound (8.7). It remains to obtain the bound for t)\\ H i( BRo y 

We deduce this bound from the interior Schauder estimates [23, Thm VI. 5] for an elliptic 
operator A: 

U(;t)\\ HH B Ro) < C (\\A^(.,t)\\H- H B Ro+l) + \m;t)\\L H B Ro+l) ) ■ (10.7) 
We use (10.4) and (10.6) with R + 1 instead of Ro in the right hand side of (10.7) and obtain 

U(;t)\\m { B Ro) < C^R, J Ro)e- at ||^|| , (i?) , * > T(R,R + 1). (10.8) 

Therefore, (10.5), (10.8) imply (8.7). □ 



11 Appendix. Fourier transform calculations 

We consider dynamics and CFs of the solutions to the system (3.2). Let F : w h- > w denote 
the FT of a tempered distribution w G S"(IR n ) (see, e.g. [6]). We also use this notation for 
vector- and matrix-valued functions. 



11.1 Dynamics in the Fourier space 

In the Fourier representation, the system (3.2) becomes Y(k,t) = Ao(k)Y(k,t), hence 



Y(k, t) = g t (k)Y (k), g t (k) = exp(Ao(k)t). 



11.1) 



Here we denote 



A (k) 



( 



1 \ 



0/ 



Qt(k) 



cos \k\t 



sin \k\t \ 
\k\ 



11.2) 



V —\k\ sin \k\t cos \k\t J 
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11.2 C variance functions in Fourier space 

Translation invariance (1.4) implies that in the sense of distributions 

E(Y (k) ® c %(k') ) = F x ^ k F y ^q (x - y) = {2n) n 5{k + k')q (k), (11.3) 

where ®c stands for tensor product of complex vectors. Now (11.1) and (11.2) give in the 
matrix notation, 

E(Y(k,t)® c Y(k',t) ) = (2n) n 5(k + k%(k)q (k)g' t (k). (11.4) 

Therefore, 

qt(x-y) := E(Y{x,t) ® c Y{y,t)) = F- l g t {k)qv{k)g' t {k). (11.5) 

References 

[1] C. Boldrighini, R.L. Dobrushin, Yu.M. Sukhov, Time asymptotics for some degenerate 
models of evolution of systems with an infinite number of particles, Technical Report, 
University of Camerino, 1980. 

[2] C. Boldrighini, A. Pellegrinotti, L. Triolo, Convergence to stationary states for infinite 
harmonic systems, Journ. Stat. Phys. 30 (1983), 123-155. 

[3] D.D. Botvich, V.A. Malyshev, Unitary equivalence of temperature dynamics for ideal and 
locally perturbed fermi-gas, Comm. Math. Phys. 91 (1983), no. 4, 301-312. 

[4] T.V. Dudnikova, A.I. Komech, E.A. Kopylova, Yu.M. Suhov, On convergence to equi- 
librium distribution, I. Klein-Gordon equation with mixing, submitted to Comm. Math. 
Phys. 

[5] T.V. Dudnikova, A.I. Komech, H. Spohn, On convergence to statistic equilibrium in 
two-temperature problem for wave equations with mixing, preprint no. 26, Max-Planck- 
Institute for Mathematics in Sciences, Leipzig, 2000 (http://www.mis.mpg.de). 

[6] Yu.V. Egorov, A.I. Komech, M.A. Shubin, Elements of the Modern Theory of Partial 
Differential Equations, Springer, Berlin, 1999. 

[7] D. Eidus, The principle of limit amplitude, Russ. Math. Surv. 24 (1969), no. 3, 97-167. 

[8] M.V. Fedoryuk, The stationary phase method and pseudodifferential operators, Russ. 
Math. Surveys 26 (1971), no.l, 65-115. 

[9] I.I. Gikhman, A.V. Skorokhod, The Theory of Stochastic Processes, Vol.1, Springer, Berlin, 
1974. 

[10] L. Hormander, The Analysis of Linear Partial Differential Operators III: Pseudo- 
Differential Operators, Springer- Verlag, 1985. 

[11] LA. Ibragimov, Yu.V. Linnik, Independent and Stationary Sequences of Random Vari- 
ables, Wolters-Noordhoff, Groningen, 1971. 



25 



[12] A.I. Komech, Stabilisation of statistics in wave and Klein-Gordon equations with mixing. 
Scattering theory for solutions of infinite energy, Rend. Sem. Mat. Fis. Milano 65 (1995), 
9-22. 

[13] P.D. Lax, R.S. Phillips, Scattering Theory, Academic Press, NY, 1967. 

[14] V.P. Mikhailov, Partial Differential Equations, Mir, Moscow, 1978. 

[15] V.V. Petrov, Limit Theorems of Probability Theory, Clarendon Press, Oxford, 1995. 

[16] N.E. Ratanov, Stabilisation of statistical solutions of second-order hyperbolic equations, 
Russian Mathematical Surveys 39 (1984), no.l, 179-180. 

[17] N.E. Ratanov, Asymptotic normality of statistical solutions of the wave equation, Mosc. 
Univ. Math. Bull. 40 (1985), no.4, 77-79. 

[18] N.E. Ratanov, Stabilisation of Statistical Solutions of Second-Order Hyperbolic Equations, 
PhD Thesis, Moscow State University, 1984. 

[19] M. Reed, B. Simon, Methods of Modern Mathematical Physics III: Scattering Theory, 
Academic Press, New York, 1979. 

[20] M.A. Rosenblatt, A central limit theorem and a strong mixing condition, Proc. Nat. Acad. 
Sci. U.S.A. 42 (1956), no.l, 43-47. 

[21] A.G. Shuhov, Yu.M. Suhov, Ergodic properties of groups of the Bogoliubov transforma- 
tions of CAR C*-algebras, Ann. Phys. 175 (1987), 231-266. 

[22] H. Spohn, J.L. Lebowitz, Stationary non equilibrium states of infinite harmonic systems, 
Comm. Math. Phys. 54 (1977), 97-120. 

[23] B.R. Vainberg, Asymptotic Methods in Equations of Mathematical Physics, Gordon and 
Breach, New York, 1989. 

[24] M.I. Vishik, A.V. Fursikov, Mathematical Problems of Statistical Hydromechanics, Kluwer 
Academic Publishers, Dordrecht, 1988. 

[25] L.R. Volevich, B.P. Paneyakh, Certain spaces of generalised functions and embedding 
theorems, Russian Math. Surveys 20 (1965), no. 5, 1-73. 



26 



